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ABSTRACT 

CO and H I velocity fields of seven nearby spiral galaxies, derived from radio-intcrferometric observa- 
tions, are decomposed into Fourier components whose radial variation is used to search for evidence of 
radial gas flows. Additional information provided by optical or near-infrared isophotes is also consid- 
ered, including the relationship between the morphological and kinematic position angles. To assist 
in interpreting the data, we present detailed modeling that demonstrates the effects of bar streaming, 
inflow, and a warp on the observed Fourier components. We find in all of the galaxies evidence for 
either elliptical streaming or a warped disk over some range in radius, with deviations from pure circu- 
lar rotation at the level of ^^20-60 km s~^. Evidence for kinematic warps is observed in several cases 
well inside i?25- No unambiguous evidence for radial inflows is seen in any of the seven galaxies, and 
we are able to place an upper limit of ~5-10 km s~^ (3-5% of the circular speed) on the magnitude 
of any radial inflow in the inner regions of NGC 4414, 5033 and 5055. We conclude that the inherent 
non-axisymmetry of spiral galaxies is the greatest limitation to the direct detection of radial inflows. 
Subject headings: galaxies: kinematics and dynamics — galaxies: spiral — galaxies: ISM 



1. INTRODUCTION 

The potential importance of radial gas flows for un- 
derstanding the evolution of galaxies has long been 
recognized. Their effect on galactic morphology has 
been discussed in the context o f evolution along 
the H ubble sequence fe.g.. iNorman. S ellwood. fc Hasan 
Il996^ and the formation of exponential stel lar disks 
ijLin Prind^[T987t IP^rgnson Cla,7k3l200l . Their 
possible role in fueling nuclear activity and star forma- 
tion in the inner disks of galaxies has been discussed by 
l^losman. Frank. & Begclman (1989) and Bhtz (1996) 
respectively. Their ability to account for the steeper than 
expected ab undance gradients in spiral galaxi es has been 
explor ed bv iLacev fc Falll lll985t). IC hamcha m fc Tavie3 
((i^gl), and'Portin ari fc Chiosil l|200a) . Finally, the sus- 
ceptibility of galactic disks to rad ial flo ws has been dis- 
cussed by iStruck-Marceiil (|199H) and iStruck fc SmithI 
l)1999D from a hydrodynamic standpoint. In general such 
flows are expected to be mildly subsonic, i.e. up to a few 
km s~^ in typical galaxies. 

Despite the strong motivations for expecting radial 
flows to occur, published measurements of radial gas 
flows in galaxy disks are limited to a few stro ngly barred 
galaxies (e.g., NGC 7479 , QuiUc n et al.ll995t NGC 1530, 
lR,egan. Vogel. fc Teubenil997ll . and even in these cases it 
is difficult to measure inflow rates unambiguously. Part 
of the problem is achieving adequate velocity resolution: 
even slow radial flows can be important on an evolution- 
ary timescale, since 1 km s^^ « 1 kpc Gyr~^. More fun- 
damental is the ambiguity in deriving the full velocity 
structure from only the line-of-sight component of mo- 
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tion. For gas in pure circular rotation, the direction of 
the velocity vector at any point in the disk is uniquely de- 
termined (assuming the orientation of the disk is known) 
and its magnitude depends only on the galactocentric 
radius. If non-circular motions exist, however, the direc- 
tion of local motion is generally not known, and recov- 
ering the true velocity structure from Doppler velocities 
alone is impossible without further assumptions. 

One approach for measuring inflow speeds is to model 
the gravitational potential based on the observed stel- 
lar distribution and then use the model to predict the 
gas inflow rate towards the center. While such modeling 
is becoming less daunting computationally, the results 
still depend quite heavily on the exact techniques em- 
ployed and on assumptions regarding the mass-to-light 
{M/L) ratio and the deprojection of surface brightness 
into a volume mass density. Determining the appropriate 
M/L ratio is especially problematic when observations of 
the stellar velocity dispersion are unavailable. An upper 
limit to M/L is provided by the "maximum disk" hy- 
pothesis, which adopts the maximum M/L ratio for the 
disk that is consistent with the rotation curve and ob- 
served gas and stellar profiles. However, the maximum 
disk hypothesis ha s been questione d, both in the context 
of the Milky Way (lMerrifielJl992t but seelSackettlll997j) 
as well as for other galaxies (|Bottemalll993fl . An alter- 
native is to determine the M/L ratio from stellar popu- 
lation synthesis mo dels. For th e strongly barred galaxy 
NGC 7479, Quille n et all (|199 5') found a consistent value 
of M/Lk ~ 1.35 Mq/ Lkq using both maximum disk 
and population synthesis techniques. By estimating the 
torque from the resulting gravitational potential, they 
inferred a radial infiow speed of 10-20 km s^^ for molec- 
ular gas along the bar. However, they only integrated 
the torque across a linear dust lane seen in CO emission; 
integrating over the entire gas distribution would likely 
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Table 1 

Properties of the Sample Galaxies 



^25^ Vq^ Nuc. CO Resolution HI Resolution 

Name R.A. (2000) Decl. (2000) (arcsec) (km s"!) Type<^ (arcsec X arcsec X km s ^) 



NGC 4321 


12:22:54.9 


15:49:21 


220 


1571 


T2 


8.85 


X 


5.62 


X 


10 


20.0 


X 


20.0 


X 


20.6 


NGC 4414 


12:26:27.1 


31:13:24 


110 


716 


T2: 


6.53 


X 


4.88 


X 


10 


17.5 


X 


15.3 


X 


10.3 


NGC 4501 


12:31:59.2 


14:25:14 


210 


2281 


S2 


7.89 


X 


6.23 


X 


20 


23.2 


X 


22.3 


X 


10.3 


NGC 4736 


12:50:53.1 


41:07:14 


340 


308 


L2 


6.86 


X 


5.02 


X 


10 


15.0 


X 


15.0 


X 


5.15 


NGC 5033 


13:13:27.5 


36:35:38 


320 


875 


SI. 5 


7.52 


X 


6.49 


X 


20 


19.9 


X 


17.0 


X 


20.6 


NGC 5055 


13:15:49.3 


42:01:45 


380 


504 


T2 


5.69 


X 


5.30 


X 


10 


12.9 


X 


12.8 


X 


10.3 


NGC 5457 


14:03:12.5 


54:20:55 


870 


241 


H 


6.88 


X 


6.36 


X 


10 


15.0 


X 


15.0 


X 


5.15 



''Semimajor axis at 25.0 mag arcsec , from RC3. 

'^Heliocentric velocity, from NED. 

•^Nuclear spectral class and type, from lHo et al J 119971) . "S" refers 
to Seyfert, "H" to H II, "L" to LINER, and "T" to transition 
(H II/LINER). Uncertain classifications are marked by a colon (:). 



result in a lower inflow speed, in better agreem ent with 
numerical simulations (e.g.. lAthanassoulalll992) . 

In this study, we will instead attempt to detect ra- 
dial inflows directly, from the atomic and molecular gas 
kinematics of nearby, moderately inclined spiral galax- 
ies. We consider three basic types of deviations from a 
circularly rotating disk: uniform (axisymmetric) inflow, 
elliptical streaming in a non-axisymmetric potential (in- 
cluding the special case of spiral streaming) , and circular 
orbits that are inclined with respect to the main disk (an 
outer warp). Since the latter two situations are not nec- 
essarily accompanied by net inflow, we aim to distinguish 
them from the pure inflow case by employing a Fourier 
analysis of the velocity field. Optical and near-infrared 
images are used as an independent indicator of possible 
bars or oval distortions. Although our models are sim- 
plistic, they have very few free parameters and hence can 
be applied to a large number of galaxies with relatively 
little fine-tuning. 

This paper is organized as follows. In fj^we summa- 
rize the observational data, most of which come from CO 
mapping as part of t he BIMA Survey o f Nearby Galaxies 
ijRegan et al.ll2001t iHelfer et al]l200a) and from previ- 
ously published VLA H I data. In ^we present velocity 
fields, tilted- ring modeling and isophote fits for the seven 
galaxies in our sample. In 21 describe the signatures 
of different types of non-circular motions, as revealed by 
a harmonic decomposition of the velocity field. In |5l 
we analyze the velocity fields for evidence of radial gas 
flows, making use of complementary information from 
the optical and near-infrared isophote fits. |H1 contains a 
discussion and summary of this work. 

2. OBSERVATIONS 

2.1. Properties of the Sample 

The sa mple of seven galaxies used in this study is the 
same as m iWong fc Blit3 ll200l hereafter Paper I). Basic 
properties of the galaxies are summarized in Tabled All 
exhibit strong CO emission at a range of radii (out to at 
least R ~ 50") and have available H I maps at a resolu- 
tion of ~20" or better. They have also been selected to 
possess fairly regular velocity fields — none are involved in 
strong tidal interactions (although NGC 4321 and 4501 
are members of the Virgo cluster) and none are classified 



as st rongly barred (SB) in RC3 l)de Vaucouleurs et alJ 
IT99l|) . Only two of the galaxies, NGC 4501 and 5033, 
contain Seyfert nuclei; th e rest sho w H II, LINER, or 
composite nuclear spectra l|Ho et al.l [T997'). At a median 
distance of 16 Mpc, adopted for the Virgo cluster us- 
ing the Ccpheid measurement of Fcrrarcsc et al. (1996), 
1" = 80 pc, and 1 kpc = 13". 

The exclusion of strongly barred galaxies may seem 
puzzling given that the gravitational torque exerted by a 
bar is considered the most effective mechanism for tun- 
neling gas into the central regions (e.g., Combes 1999). 
Many strongly barred galaxies exhibit dust lanes (usu- 
ally along the leading side of the bar) which are thought 
to be associated with shocks, and hydrodynamical sim- 
ulations indicate that these shocks can lead to strong 
radial inflows ( Athanassoula 1992). Disentangling ra- 
dial inflow from strong bar streaming, however, requires 
the type of detailed modeling for which our method is 
unsuitable, since our elliptical streaming models assume 
non- intersecting gas orbits, a condition which would not 
be satisfied where shocks occur. Hence we will focus on 
unbarred or weakly barred galaxies, which might still ex- 
hibit detectable (albeit smaller) radial flows. 

2.2. CO and HI Observations 

CO (1-0) observations were conducted with the BIMA^ 
interferometer at Hat Creek, California and the NRAO'^ 
12 m telescope at Kitt Peak, Arizona, mostly as part of 
the BIMA Survey of Nearby Galaxie s (BIMA SONG) 
ijRegan et alJ I2001I IHelfer et all 1200^ . Details of the 
observations and reduction procedures are given in Pa- 
per I. The BIMA visibilities were used to generate data 
cubes with synthesized beams of '-^6"(see Table The 
velocity channel increment was 10 km s^^ for most of 
the cubes, but 20 km s'^ for the NGC 4501 and 5033 
cubes. Although the instrumental velocity resolution was 
~4 km s^^, choosing a larger velocity channel width can 
significantly improve the signal-to-noise of the map, al- 
lowing for better deconvolution. The 12 m data were 

The Berkeley-Illinois-Maryland Association is funded in part 
by the National Science Foundation. 

^ The National Radio Astronomy Observatory is a facility of the 
National Science Foundation, operated under cooperative agree- 
ment by Associated Universities, Inc. 



Radial Gas Flows in Spiral Galaxies 



3 



folded into the BIMA cubes using the MIRIAD task IM- 
MERGE, resulting in combined datacubes that are sensi- 
tive to emission on both large and small scales. Without 
the single-dish data, the velocity field would be biased 
more strongly to the emission peaks. 

The H I observations were obtained at the NRAO Very 
Large Array (VLA). Six of the galaxies have appeared 
previously in the lit erature, in papers by lKnapen et al.l 
Ji993) (NGC 4321) . lThcmilev fc Mund^ (1997b) (NGC 
4414) IBraiml ^9&} (NGC 473 6 and 5457 ), Thcan ct ^ 
(dla) (NGC 5033), and .Thornlev fc Mundv (1997j^ 
(NGC 5055) . Reduced datacubes for these galaxies were 
kindly provided by these authors. The NGC 4501 data 
are described in detail below. All galaxies were observed 
in the C configuration, which gives a synthesized beam of 
~13" for uniform weighting and ~20" for natural weight- 
ing; three of the galaxies (NGC 4736, 5055, and 5457) 
were also observed at higher resolution in the B configu- 
ration. Data obtained in the D configuration are included 
for all galaxies except NGC 4501 and 5033, so that the re- 
sulting maps are sensitive to even very extended (~15') 
emission. The datacubes were imported into MIRIAD 
and transformed to the spatial and velocity coordinate 
frames of the BIMA cubes. 

2.3. NGC 4501 Data 

The single-dish CO map for NGC 4501 used in Paper 
I was only 2' x 2' in size and had an angular resolution 
of 55", making it less than ideal for combination with 
the BIMA data. For this study we have instead used a 
CO map of NGC 4501 taken during several sessions in 
1992 and 1993 with the IRAM 30 m telescope on Pico 
Veleta, Spain. The map consisted of spectra taken at 249 
positions separated by 12", covering a field roughly 3'x 
2' elongated along the major axis. At A = 2.6 mm the 
half-power beamwidth (HPBW) of the 30 m telescope is 
23". The noise level in the resulting map was 27 mK per 
10.4 km s~^ channel. The data were converted from a 
brightness temperature to flux density scale assuming a 
telescope gain of 6.5 Jy (K[T^])-i. The resuhing CO flux 
within the primary beam function of the BIMA mosaic 
is 2200 Jy km s~^, comparable to the 2300 Jy km s~^ 
measured in the 12 m map. We used the IRAM map to 
regenerate a combined (single-dish + interferometer) CO 
datacube for this galaxy using the IMMERGE technique 
as described in Paper I. In performing the combination, 
the flux of the BIMA map had to be rescaled by a factor 
of 0.8 to produce consistency with the IRAM map. This 
scaling factor is similar to factors needed to reconcile the 
BIMA and 12-m flux scales for the other galaxies (Paper 
I), and is consistent with the calibration uncertainties 
typical of millimeter-wave telescopes. 

As the H I data for NGC 4501 used in this study (as 
well as Paper I) have not previously been published, we 
also briefly describe the observations and data reduction 
here. The observations were made on 1991 January 20 as 
part of a program (AG 318) designed by Gunn, Knapp, 
van Gorkom, Athanassoula and Bosma to determine ro- 
tation properties of spiral galaxies. The instrumental pa- 
rameters are given in Table[21 Standard VLA calibration 
procedures were applied in AIPS. The continuum was 
subtracted by making a linear visibility fit for a range 
of line free channels at either side of the band. The line 
free channels were determined by visual inspection of the 



Table 2 

Parameters For NGC 4501 H I Observations 



R.A. of phase center (1950.0) 12''29™2r;0 

Decl. of phase center (1950.0) 14°41'43'.'0 

Central velocity (km s^^, heliocentric) 2280 

Velocity range (km s^^) 630 

Time on source (hr) 5.8 

Bandwidth (MHz) 3.125 

Number of channels 63 

Channel separation (km s~^) 10.3 

Synthesized beam width 23'.'2 X 22'.'9 
Noise level (Icr) 0.55 mjy beam~^ 



data. The line images were made with the task IMAGR, 
using a "robust" weighting scheme intermediate between 
uniform and natural weighting, as well as a tapering of 
the longer baselines. This resulted in a synthesized beam 
of 23"2 X 22"9 (FWHM), close to the beam size of the 
30-m CO data. 

2.4. NIR and Optical Images 

Near-infrared (NIR) images were obtained from the 
Two Micron AU Sky Survey (2MASS) in the (2.2 
/xm) filter for all galaxies e xcept NG C ^41 4 and 5055, 
for which K' images from iThornlevI (l996j) were used 
instead . Astrometry and photometry from 2MASS and 
iThornl ev (1996:1) were adopted; surface photometry from 
the latter was found to be consistent with 2MASS for 
NGC 4414 to within <0.5 mag arcsec"^. Astrometric co- 
ordinates should be accurate to less than 1", as indicated 
by a comparison with Digitized Sky Survey (DSS) images 
using foreground stars. In some cases (NGC 4501, 4736, 
5033, 5457), the galaxy feU near the edge of a 2MASS 
Atlas Image, and several Atlas images were mosaiced to- 
gether using the IRAF task COMBINE, after subtracting 
a constant background level from each. 

We also obtained /-band images for all galaxies except 
NGC 4736. For NGC 4321, 4414, 4501, 5033, and 5055 
these images come from the on-line catalog of Frci et al] 
l)1996tl . For NGC 5457, we used an / image taken with 
the 1.5 m Palomar telescope in 2000 April, taken by M. 
Regan for the BIMA SONG project. For NGC 4736, 
we used instead the A=665 nm continuum image corre- 
sponding to the Ha image of Martin & Kcnnicutt (2001), 
approximating an i?-band image. As most of these im- 
ages were not photometrically calibrated, and some are 
saturated near the galaxy's nucleus, they are only used 
her e to sup plement the if-band data. Astrometry for 
the iFrei et al. (1996) images proved difficult, since fore- 
ground stars have been removed, and was accomplished 
by matching the location of giant H II regions with a DSS 
image (NGC 4321, 5033, 5055) or by a cross-correlation 
technique (NGC 4414, 4501). We expect the resulting 
coordinates to be accurate to within 2". 

3. DATA ANALYSIS 

3.1. Velocity Fields 

We derived CO and H I velocity fields from the corre- 
sponding datacubes by fitting a Gaussian to the spec- 
trum at each pixel using a customized version of the 
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Fig. 1. — Velocity fields for all seven galaxies derived from Gaussian fits as described in the text. For each galaxy, the left panel shows 
the H I velocity field, the center panel shows the full-resolution CO velocity field, and the right panel shows a contoured / or ii'-band image 
(note the presence of some foreground stars). Units of the axes are in arcseconds. The velocity fields have been masked using a rotation 
model as described in the text. Isovelocity contours are spaced by 20 km s~^. The synthesized beam is shown on the lower left, and the 
CO field of view (where the map sensitivity is at least half of its peak value) is given by a dotted contour. 



Radial Gas Flows in Spiral Galaxies 



5 





6 



MIRIAD task GAUFIT. Gaussian fits were required to 
have an ampHtude of >1.5cr, where a is the r.m.s. noise 
in a channel map. Although this criterion allows for the 
detection of broad lines which are relatively weak in a sin- 
gle velocity channel, it has the drawback of also allowing 
for a large number of spurious fits, as evidenced by veloc- 
ities which are spatially uncorrelated from pixel to pixel, 
producing a "noisy" appearance. In contrast, the true 
velocity field should be dominated by galactic rotation 
and smooth on scales smaller than a synthesized beam. 
Additional thresholds were therefore imposed based on 
the integrated flux of the Gaussian and the uncertainty 
in its mean velocity. Although the exact threshold de- 
pended somewhat on the signal-to-noise and flux level 
in the datacube, a Gaussian fit in a naturally weighted 
CO datacube {a ~ 50 mJy bm^^) would typically be 
required to have a velocity-integrated intensity of > 3 
Jy bm~^ km s~^ (a 3a detection over 20 km s~^) and a 
velocity uncertainty of < 20 km s~^ to be included. The 
assumption of Gaussian profiles, of course, is only an ap- 
proximation, but in most cases the signal-to-noise could 
not justify the use of an asymmetric profile, and other 
methods of deriving the velocity field (such as the first 
moment) proved less robust and did not provide mean- 
ingful error estimates. 

Figure shows the CO and H I velocity fields at full 
resolution and a contoured / or K-hand image on the 
same scale as the CO velocity field for comparison. The 
BIMA field of view (where the sensitivity drops to half 
of its peak value) is shown as a dotted contour. The 
velocity fields shown have been masked further using a 
rotation model derived from the adopted rotation curve. 
This masking eliminates points which deviate by >50 
km s^^ from circular rotation, as described in H3.2I As 
noted in Paper I, the CO emission in all galaxies except 
NGC 4736 and 5033 appears to extend to the edge of the 
observed field, suggesting that our observations sample 
only the inner part of the molecular disk. 

Note that while there is some subjectivity in determin- 
ing the appropriate thresholds for accepting the Gauss- 
sian fits, our velocity fields do not automatically reject 
anomalous velocity gas that is inconsistent with circular 
rotation, until the final masking using a rotation model 
is performed. Indeed an anomalous velocity component 
is seen in the southeastern corner of the H I disk in NGC 
4501. No other clear examples of anomalous velocity gas 
were seen. 

3.2. Kinematic Disk Parameters 

In preparation for the more detailed kinematic anal- 
ysis needed to search for radial fiows, we estimated the 
basic disk orientation parameters (mean velocity, rota- 
tion center, position angle, and inclination) as well as 
the rotation curve by fitting tilted-ring models to the 
CO and H I velocity fields. Four velocity fields were 
analyzed for each galaxy: the CO and H I at full resolu- 
tion, the CO data smoothed to the resolution of the H I, 
and the H I data smoothed to 55" resolution. We used 
the stand a rd lea st-squares fitting technique developed by 
iBeeemanl 1)1987 !). as implemented in the ROTCUR task 
within the NEMO software package (Teubcn 1995). The 
galactic disk is subdivided into rings, each of which is de- 
scribed by 6 parameters: the center coordinates (xo,?/o)i 
position angle (Fq), inclination (i), offset velocity (T4ys), 



Table 3 

Galaxy Parameters From Kinematic Fits 



Galaxy 




V b 

► sys 


To" 


i 




(") 


(km s-l) 


(°) 


n 


NGC 4321 


(0,0) 


1570 ± 3 


153 ± 3 


34 ± 5 


NGC 4414 


(2,0) 


720 ± 5 


159 ± 2 


55 ± 2 


NGC 4501 


(0,0) 


2270 ± 1 


141 ± 1 


64 ±2 


NGC 4736 


(1,1) 


317 ± 4 


295 ± 10 


32 ±8=^ 


NGC 5033 


(-1,0) 


880 ±8 


353 ± 2 


68 ±2 


NGC 5055 


(0,0) 


506 ±3 


98 ± 1 


63 ± 1 


NGC 5457 


(0,0) 


257 ± 1 


42 ± 2 


21 ± 3 



''Position of kinematic center with respect to optical center, in 
directions of RA and DEC. 
''Systemic LSR velocity. 

"^Position angle of receding side of line of nodes, measured E from 
N. 

'^Inclinatio n of 35° was a dopted based on the photometric study 
of iMoUenhoff et al.l 11993) . 



and circular velocity (wc). Starting with initial estimates 
for the fitting parameters, the parameters are adjusted 
iteratively for each ring until convergence is achieved. 

In general it is not possible to fit all 6 parameters si- 
multaneously, due to kinematic disturbances and limited 
signal-to-noise in the velocity field, compounded by cor- 
relations between fit parameters. We therefore adopted 
the following procedure for determining the fit parame- 
ters: 

1. Initial estimates for i and Fq were taken from the 
isophotal fits ( H3.3|l . and the optical or radio contin- 
uum center from the NASA Extragalactic Database 
(NED) was taken as {xqjI/q). Initial estimates for 
V^ys and Vc were taken from a position-velocity cut 
through the center at a position angle of Fq, or 
from previously published studies. 

2. An improved estimate of Fq was determined by fix- 
ing (xq, yo), Vsys, and i, and allowing the remaining 
parameters (Fq and Vc) to vary. 

3. An improved estimate of V^ys was determined by 
fixing (xo,?/o), To, and i. 

4. An improved estimate of (xo,yo) was determined 
by fixing Vsys, Fq, and i. 

5. An attempt was made to fit for i while (xq, 2/0)7 
V^ys, and Fq were held fixed, although this could 
only be accomplished for some of the rings. 

At each step, an inner radius for the fit was chosen to 
exclude rings within 2 beamwidths of the center, which 
were strongly affected by finite resolution effects (beam 
smearing). An outer radius was chosen to ensure that the 
ring was well-sampled in azimuth and did not show signs 
of a strong warp (a systematic drift in Fq with radius). 
Since the program fits each ring separately, a weighted 
average of the values for all rings was used as the im- 
proved estimate for the next step. The full-resolution 
H I maps were generally most useful for constraining the 
fit parameters, since they contain a large number of in- 
dependent resolution elements in a region where Vc is 




roughly constant, but the CO maps provide better con- 
straints on the kinematic center since they have higher 
resolution (and often signal-to-noise) in the central re- 
gion. 

Once an initial rotation curve and set of orientation 
parameters had been determined, model velocity fields 
were generated (assuming a flat rotation curve beyond 
the outer edge of the fitted disk) and compared to the ob- 
served velocity fields. Pixels in the velocity fields which 
were discrepant by more than ~50 km s~^ from the cor- 
responding model were flagged. This threshold was in- 
creased when necessary so that only points that were 
clearly noise were rejected: care was taken to ensure that 
the threshold did not exclude gas that was continuous 
in velocity with gas in the disk. Although this process 
introduces some noise bias (noise pixels with velocities 
within 50 km s~^ of the model are not eliminated), it pro- 
vides an excellent compromise between sensitivity (which 
would be degraded by too strict a Gaussian fit thresh- 
old) and noise rejection, as evidenced by the generally 
smooth velocity contours in Fig. ^ 

The cleaned velocity fields were then subjected to the 
ROTCUR analysis again to derive "final" values for the 
fit parameters, which are summarized in Table |3 The 
uncertainties given represent the range in values for the 
various velocity fields analyzed, and are thus generally 
larger than the formal least-squares errors for a single ve- 
locity field. In ^4. 51 we discuss how errors in the adopted 
disk parameters would affect the resulting harmonic de- 
composition. 

Note that in most of the galaxies the kinematic center 
is displaced < 1" from the center defined by the K-h&nd 
or radio continuum nucleus. Only NGC 4414 and 4736 
show a displacement of ~1'.'5, which is of marginal signif- 
icance given the uncertainties in determining the optical 
and kinematic centers. For the subsequent analyses pre- 
sented here, we have adopted the kinematic rather than 
optical center. 



3.3. Isophote Fits 

In order to better constrain the influence of spiral 
structure and bars on the observed kinematics, we per- 
formed ellipse fitting to the optical and NIR isophotes 
using the IRAF task ELLIPSE. The fitting program out- 
puts the center {x, y) coordinates, semimajor axis a, el- 
lipticity e = 1 — h/ a, and position angle of each fit. The 
ellipse centers were fixed during the fitting except in the 
case of NGC 4321, for which a satisfactory fit could not 
be obtained without allowing the center to wander by up 
to ^10". However, the difficulty in fitting the isophotes 
does not reflect an uncertainty in the center position: the 
nucleus of NGC 4321 is well-determined photometrically 
in the if-band image (to within <2") and the kinematic 
center is also well-determined (to within '-^l") by the CO 
velocity field. The results of the isophotal analysis are 
discussed in tl5.2l 

4. SIGNATURES OF NON-CIRCULAR MOTIONS IN DISKS 

In this section we present some general results for 
non-circular motions in disks. We adopt a right- 
handed xyz coordinate system with the main disk of 
the galaxy confined to the x y plane. As discussed by 
iFranx. van Gorkom. fc de Ze euw ( 19!^, the direction of 
our sight-line is defined by two angles, 6'obs and i, where 
^obs is measured in the xy plane from the x-axis and i is 
measured from the z-axis [see Figure l^^a)]. Note that i 
is restricted to the range < i < 7r/2. Figure |2fb) then 
shows the appearance of the galaxy as viewed in the sky. 
Here F is the position angle measured in the sky plane, 
from north to east, and Fnodo is the position angle of 
the line of nodes (the intersection of the galaxy and sky 
planes). Although in principle we could define the x-axis 
so that 0obs=O, we choose to keep it independent of 6'obs 
to more easily handle the case of a bar. 

Adopting the convention that positive line-of-sight 
("radial") velocities correspond to recession, we there- 
fore have: 

V\os = Vsys + Vc cos{e ~ 9ohs - 7r/2) sini (1) 
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for an axisymmetric model with rotation curve Vc{R). 
Here Vsys is the mean (systemic) velocity of the galaxy. 
For the most general case of a two-dimensional velocity 
field, 

Vios{x,y) = Vsys+vg{x,y)cosil;smi+vji{x,y)sm'4jsini , 

where (vR^ve) is the velocity vector in the disk in po- 
lar coordinates and ^ = — 9ohs T '''/2 is the azimuthal 
angle in the galaxy plane, measured from the receding 
side of the line of nodes (the — sign is taken for counter- 
clockwise rotation and the -I- sign for clockwise). With 
these definitions vg is always positive, and vr > in- 
dicates outflow for counterclockwise rotation and inflow 
for clockwise rotation. 

In each case below, we seek to express the model ve- 
locity field Vios in terms of harmonic coefficients Cj{R) 
and Sj{R), where 



First we consider coUisionless models, which are 
strictly applicable only to stellar orbits. Using the po- 
tential in Eq. |31 and the equation of motion in a frame 
rotat ing at angular velocity flp (Binncv & Trcmain^ 
one can solve for the allowed closed, pe- 
riodic orbits using first-order epicyclic theory (this 
assumes the bar perturbation i s rel atively weak). 
iSchoenmakers. Franx. fc de Zeeuwl l)1997|) showed that a 
potential distortion of order m introduces changes in the 
(m — l)st and (to + l)st harmonic coefficients. For the 
m—2 case, their expressions reduce to: 



sys 



(6) 



ci = <j 1 - - [(2 - 3ujp)a2 + (1 - Wp)&2] cos 26*0^ 



si = <j - [(2 - 3wp)a2 + (1 - ujp)b2] sin 29ohe 



los 



Co 



[cj cos{jip) + Sj sin(j?/>)] 



(3) 



C3 = w* |--[(wp - 2)a2 + (1 - LUp)b2] cos26'obs 
S3 = w* I^K'^p - 2)a2 + (1 - cjp)&2] sin26'obs| 



rip/Qo — R{Qp/vc), and the 



where v■^, = VcSini, uj 

ellipticity coefficients 02 and 62 are given by: 

1 



a2 



(1 - iuj 

2 - 3lu, 



p , 2L0l)il~iUp) 



Epot 



(7) 



(l-4Lc;p + 2cj2)(i_^^)2 



tpot 



For the simple models we consider, it is only necessary to 
expand out to n—3. For the case of pure circular rotation, 
co=ysys and ci=WcSini, with all other coefficients 0. 

4.1. Axisymmetric Radial Flow 

The simplest example of non-circular motion is one in 
which vn and vg are functions of radius R only: 

V\os = Vsys + vg{R) COS V' sin i -t- Vfj{R) sinipsini (4) 

This corresponds to circular motion superposed on ax- 
isymmetric inflow or outflow. The harmonic coefhcients The coefficients si and S3 for this model are plotted with 
are simply co=14ys, ci=?;esini, and si=vnsmi. Of a dashed line in Figure Ela) as a function oi cUp, which 

is proportional to the radius R. Note that between the 
ILR (inner Lindblad resonance, at iVp — 1 — \/2/2) and 
OLR (outer Lindblad resonance, at = 1 + \/2/2), the 
si and S3 terms have opposite sign. 

Alth ough gas is not expecte d to follow coUisionless 
orbits. iGerhard fc V ictrT (19^^ and lFranx et alJ l|1994^ 
have argued that a coUisionless treatment should be ap- 
plicable to gaseous orbits if the potential is non-rotating, 
since gas will then be able to settle into stable closed or- 
bits. The special case of a stationary pote ntial {ujp=0) 
was solved previously by iFranx et alJ l)1994j) : 



course, for the inflow case a continuity problem arises 
unless some sink for the inflowing gas exists: this might 
be star formation, a central black hole, or an outflowing 
wind ejected out of the plane. 

4.2. Elliptical Streaming 

A straightforward alternative to simple radial flows is 
the case in which gas follows elliptic closed orbits, as 
might occur in a bar-like potential. For a general distor- 
tion of harmonic number m rotating at a fixed pattern 
speed rip, the potential is given by: 

$(i?, 9) = <^>o{R) + 4>,„(i?) cos[m(6' - Qpt)] (5) 

where $0 is the axisymmetric part of the potential. For 
simplicity we take to be a spherical logarithmic po- 
tential, giving a flat rotation curve with circular speed 

Vc- 

^o{R)=vllnR. 

We assume a uniform to=2 perturbation to the potential 
given by: 

1 



ci = ( 1 - epot cos 26'obs ) 
si =w*epot sin 26'obs 
C3 = S3 = 



(8) 



so that the ellipticity of the potential at all radii is 
Epot- The sign of <I>2 is taken to be negative in order 
that the equipotentials are elongated along the x-axis 
{9=0} at t=0. Potentials of this type have been ana- 
lyzed in a number of pr evious studies (e.g., Franx et al 
1 1994 iKuiiken fc Tremaina il994 iRix fc Zaritskvi il99 ' 

i,Tog|i2nnnD. 



Here no resonances exist and all orbits are aligned along 
the y-axis, perpendicular to the bar. The fact that the 
C3 and S3 terms vanish has an important implication: 
for a flat rotation curve, the velocity field of a station- 
ary bar potential resembles axisymmetric inflow or out- 
flow. Thus kinematics alone cannot distinguish these two 
cases; in i |4.6l we will discuss how photometric measure- 
ments might break this degeneracy. Since measurement 
of ci only provides the value of (1 — epot cos 29ohs), and 
cannot be determined independently, it is impossible 
to disentangle epot from the viewing angle 6'obs- How- 
ever, under the assumption that ci « v^:, a measure- 
ment of the si term provides an estimate of the quantity 
epot sin20obs, and hence a lower limit on epot- 
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Fig. 3. — (a) Harmonic coefficients si and S3 normalized by Vt, = ricsini, plotted as a function of ijp oc _R for a coUisionless model 
(dashed line) and a dissipative model {solid line) with a flat rotation curve and a rotating, globally elongated potential. The locations of 
the ILR, CR, and OLR are marked as dotted lines, (b) Same but for a two-armed spiral model with a fiat rotation curve, pitch angle of 
10°, and velocity perturbation amplitude Vsp = 0.2vc- 



For the more general case of gas flow in a rotating po- 
tential, as in a barred galaxy, it is possible to relax the 
coUisionless assumption by including a damping term in 
the equation of motion to simulate the dissipative ef- 
fect of gas viscosity. This can be done as described in 
Appendix ^ The corresponding harmonic coefficients 
Si and S3 are plotted with a solid line in Figure |3J a). 
Their behavior is qualitatively similar to the coUisionless 
model, aside from a smoother variation near the reso- 
nances. For the non-rotating case, the dissipative and 
coUisionless models yield the same coefficients. 

To determine the effect of using a more realistic po- 
tential, we a lso ran models using the analytic potential 
employed bv iWadal ljl994j) . which gives a rotation curve 
that can be parametrized as Vc oc R{R^ + a^)"^/**. Al- 
though the results differ in detail from the flat rotation 
curve case, the si and S3 coefficients still have opposite 
signs except near major resonances. The region where 
the rotation curve rises is associated with a second in- 
ner Lindblad resonance (the inner ILR) across which the 
coefficients both change sign. 

4.3. Spiral Arm Streaming 

A two-armed spiral density wave can be thought of as 
an m=2 perturbation with a phase shift that varies with 
radius. One would therefore expect that the velocity per- 
turbations due to spiral arms would be in rough qualita- 
tive agreement with the results given in i j4.2l For a more 
de tailed analysis we adopt the linearized equations given 
bv lCanzian fc Alien! l(T997.) for the velocity perturbations 
due to a spiral potenti al, which are based on th e origi- 
nal treatment given bv lLin. Yuan, fc Shd llT969l) . Sub- 
stituting into their equations, we find for a two-armed 
logarithmic spiral the harmonic coefficients 

Cl - Vc _ Vsp /_K 



smi 



2 V20 



ly sin (26'sp - x) 



(9) 



^ = cos(2e.p-x) 
smi 2 \2it / 

^ = + Sin(203p + X) 

Sim 2 \2i2 / 



■S3 

sin I 



Vsp f K 

2 \2n 



cos (26lsp + x) 



Here 9sp is the spiral phase, defined at a fiducial radius 
Ro by 

ln{R/Ro) 

fsp — , 

tanx 

where x is the pitch angle, the angle an arm makes with 
the tangent to a circle. The dimensionless frequency v 
and epicyclic frequency k are defined as 



in^ + R- 



dR 



and 



is a velocity amplitude that depends on the 



strength of the spiral perturbation. 

The si and S3 coefficients for a flat rotation curve are 
plotted in FigureO^b) as a function of ojp oc R. A charac- 
teristic of the spiral wave is the sinusoidal variation in the 
coefficients with radius — since the spiral has no preferred 
orientation, the observer in effect samples different orien- 
tations of the perturbed velocity field at different radii. 

In a realistic treatment including the nonlinear ef- 
fect of shocks, the spiral velocity perturbations will not 
be purely sinusoidal (e.g., Roberts & Stewart 1987), the 
peaks in Figure |2Ib) will be narrower, and additional 
Fourier terms will occur. However, two basic properties 
of the model are worth noting since they are likely to 
be quite general. First, there is the change in dominance 
from t he si term to th e S3 term that occurs at corotation, 
which ICanzianI ljl993fl showed must occur when crossing 
a corotation resonance. The same phenomenon occurs in 
the bar model, although to a lesser degree [Figure |2Ia)]. 
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Second, the signs of si and S3 are the same near corota- 
tion (where v = 0), but tend to take on opposite signs 
as one moves away from the corotation resonance (CR). 
For the bar models, on the other hand, the two coefFi- 
cie nts have opposite signs at nearly all radii. As shown 
bv lSchoenmakers et al.l l|1997!) . this difference can be at- 
tributed to additional contributions to si and S3 that 
result from the radial variation in the phase of the po- 
tential. 

4.4. Warped Disk 

Although not directly related to radial gas flows, warps 
are a common feature in the outer H I disks of galaxies, 
and are well-known to have a significant effect on the 
observed velocity field fe.g.. iBosmal 11973 ). A warp is 
characterized by gas moving in circular orbits, but with 
the angular momentum vectors of the orbits not aligned 
as they would be in the case of a flat disk. Most galaxies 
are not strongly warped within their optic al disks (for a 
review of observed properties of warps see IBriggg| ll99ClD . 
so the plane of the inner disk defines a natural coordinate 
system for the galaxy. 

The kinematics of warped galaxies can usu- 

ally be well-modeled wi th a tilted-ring fit 

l|Rogstad. Lockhart. fc \\^gh^ ITqtI iBosmal ll98H) . 
in which the inclination i and the major axis position 
angle Fmaj are allowed to vary with radius. In contrast 
to the cases of radial flows or elliptical streaming, 
the kinematic axes Tmaj and r,nin (see H4.6(l remain 
perpendicular, even as they drift with radius ijBosmal 
^81), since the gas flow is still circular. A harmonic 
decomposition of the velocity field using fixed disk 
parameters will display the characteristic signatures due 
to incorrect values of the position angle and inclination 
discussed below f ti4.5|) . The coefficients ci, si, C3, and 
S3 will all change with radius, with the ci and Si terms 
being most affected. Again this contrasts with the case 
of elliptical or spiral streaming, where changes in the 
C3 and S3 terms will dominate beyond the corotation 
radius. 

4.5. Errors in the Parameters and Beam Smearing 

We have assumed up to this point that the disk param- 
eters used in making the harmonic expansion (namely, 
the kinematic center, inclination, and position angle) are 
correct. It is clear that this will not generally be the 
case, since the disk parameters will likely be derived from 
the v elocity field under the assu mption of circular rota- 
tion. ISchoenmakers et al] l)1997|) analyzed the changes in 
the harmonic coefficients that result from using incorrect 
disk parameters, and found that the Fourier components 
due to non-circular motions will mix with those due to 
incorrect disk parameters. To first order, and assum- 
ing the errors in the position of the kinematic center are 
small compared to the radii of interest (satisfied for our 
galaxies, see Table IS)), the following rules apply: 

1. An error in the position angle Fq leads to offsets in 
the si and S3 terms proportional to = VcSini. 
The offset in si is generally larger than that in S3. 

2. An error in the inclination i leads to offsets in the 
ci and C3 terms which are proportional to u». Be- 
cause the ci term is large, the effect will only be 
detectable in the C3 term. 



3. An error in the kinematic center leads to offsets in 
the Co, C2, and S2 terms proportional to v^/ R. 

What is the effect of finite resolution (beam smearing) 
on the coefficients? Smoothing a datacube to lower res- 
olution affects primarily the ci term (since the rotation 
curve rise becomes more gradual) and the C3 term (since 
the isovels near the galaxy center become more parallel, 
leading to a lower fitted inclination). On the other hand, 
the si and S3 terms, which are sensitive to the kine- 
matic position angle, should be relatively unaffected by 
beam smearing outside of the central few resolution ele- 
ments, since smoothing an image should not lead to any 
net rotation of the velocity contours. (Of course, any 
small-scale variations in these terms will be smoothed 
out, but there should be no systematic offsets.) We have 
confirmed these effects using simulations of axisymmetric 
velocity fields convolved wit h both symm etric and asym- 
metric beams (for details see lWongl2000|) . We note, how- 
ever, that our simulations have assumed an axisymmetric 
filled disk; if the gas distribution is very inhomogeneous, 
or confined to a small range in azimuth, it is conceivable 
that beam smearing could affect the measured si and S3 
terms as well. 

4.6. Inclusion of Isophotal Data 

Given that it is nearly impossible to kinematically dis- 
tinguish a radial inflow model from an elliptical stream- 
ing model when ujp = Qp/D,Q <C 1, i.e. when the potential 
is rotating slowly compared to the local circular speed, 
we must consider whether additional constraints can be 
provided by optical or near-infrared surface photometry. 
It is worth noting at the outset that since the isophotes 
(contours of constant surface brightness) of an axisym- 
metric disk are concentric aligned ellipses, any change 
in the isophotes over a region where the kinematics sug- 
gest pure radial flows would raise the possibility that 
streaming motions are occurring. A clearer prediction 
can only be made if one restricts consideration to an el- 
liptical streaming model in which the gas orbits are elon- 
gated in the same direction as the isophotes. Then, as 
shown below, the morphological and kinematic position 
angles will deviate in opposite directions with respect to 
the line of nodes. 

First, however, we need to identify conditions under 
which the elongation of the gas orbits will be parallel 
to the elongation of the isophotes. We consider three 
possibilities, two of which would satisfy this requirement: 

1. For a stationary elongated potential, the closed 
loop orbits are perpendicular to the elongation of 
the potential, so if the stars follow loop-type or- 
bits like the gas, the isophotes will be aligned with 
the gas orbits. In this case, however, the stars ob- 
viously cannot support the potential — in fact they 
will tend to counteract its asymmetry (Jog 2000). 
Rather, this case would be appropriate for a disk 
galaxy embedded in a massive triaxial halo. 

2. Most of the stars could follow box- type orbits 
that are aligned with the potential, so that the 
isophotes are perpendicular to the gas orbits. 
This would allow a self-consistent solution between 
the stellar density and the potential it generates 
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Fig. 4. — (a) Offset of morphological position angle from the line of nodes, as a function of Sobsi for isophotes of ellipticity 0.1. (b) 
Offset of kinematic minor axis position angle from the line of nodes, for an elongated potential with epot=0.1. In both panels, a range of 
inclinations i is shown. The numbers in parentheses in the legend give the maximum value (in degrees) for the corresponding curve. 



I)de Zeetiw. Hunter, fc SchwarzschildllQSTT. and no 

triaxial halo would be required. 

3. The potential of a rotating bar will be aligned with 
the stable loop orbits between ILR and CR (these 
are often referred to as xi orbits), as shown in Fig- 
ure IA13I in the Appendix, so again the isophotes 
will be parallel to the gas orbits as in Case 1. 

Thus, by assuming that the gas orbits are elongated par- 
allel to the isophotes, we are assuming that Case 1 or 3 
applies. Note that Case 1, which implies a strongly sub- 
maximal disk, conflicts with most observational studies 
of high surface brightness galaxies (jBell fc de J ong 2001^ 
IBottemalll993H . On the other hand, Case 3 is not ruled 
out by such considerations and should be readily appar- 
ent in the isophotes, since bars are strong between the 
ILR and CR and weaker outside this region. We assume 
throughout our discussion that the isophotes trace a thin 
disk, so we can neglect the complications of projecting 
an inherently 3-dimensional structure onto the plane of 
the sky. 

If an orbit (or isophote) has an intrinsic axis ratio 
Q = b/a and position angle ■00 with respect to the line 
of nodes, how will it appear when projected on the sky 
at some inclination angle i7 The appropriate geometri- 
cal transformations can be found inp^eubcn (1991), who 
finds that the new axis ratio q and position angle 70 (still 
measured with respect to the line of nodes) are given by 



{A cos^ i + C) cos 2V'o + {A cos^ i - C) 
{A cos2 i + C) cos 2V'o - {A cos^ i - G) 
2B cos i 

tan 270 



A cos^ i - 

where the coefficients are defined as 



C 



A = Q^ cos^ Vo + 



sin tpQ 



(10) 

(11) 

(12) 



B = 



_(i-Q') 



sin cos -00 



C = Q2 sin^V^o 



COS^ 00 



Note that we have c orrected some t ypesetting errors in 
the original paper bv lTeubenlllT99l . 

Using these equations and the sky coordinate system 
described by Figure El we can calculate the offset 70 be- 
tween the apparent or morphological position angle of 
an elliptical isophote (Fmorph) and the position angle 
of the true line of nodes (Fnodc)- For the model dis- 
cussed in ij4.2l with equipotentials elongated along the 
X-axis, we assume that the isophotes are elongated along 
the y-axis. The resulting offset angle is shown in Fig- 
ure^a) as a function of the viewing angle Qohs (measured 
from the -|-x-axis), for an ellipticity of 0.1 ((5=0.9) and a 
range of inclination angles. Note that for low inclinations 
(cosi > Q), the offset angle can reach very large values 
if the line of sight is close to the long axis of the ellipse. 

For comparison, the offset between the kinematic mi- 
nor axis (the locus in the sky plane where Vios — 



sysy 



and the line of nodes, for an orbit in an elongated po- 
tential, is plotted in Figure^b) for the same set of incli- 
nation angles (the kinematic major axis shows a similar 
but weaker effect). This angle c an be expressed in terms 
of the harmonic coefhcients as ijFranx et al.lll993) : 



node — 2 



Si - S3 1 



c\ — 3c3 cos i 



(13) 



For all non-degenerate observing angles, the offset in 
rmin is in the opposite direction from the offset in Fmorph- 
This is a robust result, depending purely on geometry, as 
can be seen by drawing on a sheet of paper an ellipse that 
is centered on, but has some angle to, the vertical and 
horizontal axes. The quadrants through which the appar- 
ent minor axis passes will be different from the quadrants 
where the maximum y-value is achieved, which is where 
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Fig. 5. — Comparison of si and S3 coefficients for (a) inflow [open squares) and warp {solid circles) models, (b) bar streaming model 
including dissipation, and (c) spiral streaming model with a pitch angle of 10°. The dashed line in panel (a) represents the "warp line" 
where points should fall if there is an error in the P.A., for the assumed inclination of 55°. Long-dashed and short-dashed lines in panels 
(b) and (c) represent fits to points inside and outside of corotation respectively (shown with different plot symbols). 



the kinematic minor axis occurs if the x-axis is the hne of 
nodes. Note also that a larger incUnation reduces r,„orph 
but increases Fmin, so bars in highly inclined galaxies are 
best identified kinematically rather than morphologically. 

5. RESULTS 

Before returning to the analysis of the seven galaxies 
in our sample, let us summarize our basic approach to 
diagnosing various types of non-circular motions. 

Based on the results of 21 the effects of elliptical or 
spiral streaming are most easily distinguished from ra- 
dial inflow by examining the first and third order sine 
coefficients (si and S3). Variations in these coefficients 
as a function of radius can be placed into three broad cat- 
egories, corresponding roughly to warp, elliptical stream- 
ing, and inflow models: 

1. The si and S3 terms are correlated, with dss/dsi > 
0. 

2. The Si and S3 terms are anti-correlated, with 
ds^/dsi < 0. 

3. The si term is significant in some region but the 
S3 term is negligible (S3 ~ 0). 

We reiterate that these criteria do not always permit a 
clear distinction between physical models. As noted in 
a stationary bar potential (usually in the second 
category) can mimic the kinematic signature of radial in- 
flow (third category) inward of the ILR. A spiral stream- 
ing model (essentially equivalent to elliptical streaming, 
but with the orbits misaligned) can produce all three 
types of behavior, depending on location relative to the 
CR, but will generally correspond to case 2. 

Examples of how four different kinematic models ap- 
pear in the S1-S3 plane are given in FigureEl For a pure 
warp model, the (si, S3) points are predicted to lie along 
a line we refer to as the "warp line," defined by 



where q = cos?. This is the predicted relation between 
si and S3 for an incorrect disk position angle f ij4.5l) . For 
elliptical streaming, the results depend somewhat on the 



choice of model parameters, but in general one finds a 
shallow slope inside of the corotation radius (Rcr) a-nd 
a steep slope outside. This is a reflection of the switch in 
dominance from the si to the S3 term as one crosses Rqr 
(Canzian Allen 1997i^- and also occurs for the spiral 
streaming model. As is clear from Figure the great- 
est ambiguity arises when trying to distinguish inflow 
from streaming motions in the region R <^ Rcr- For 
a large sample of galaxies viewed from random orienta- 
tions, streaming motions should produce apparent out- 
flow in roughly half of the galaxies, but our sample is too 
small for such statistical arguments to be useful, 

5,1, Harmonic Analysis of Velocity Fields 

Using the GIPSY task RESWRI l)Schoenmak ers'199g^. 
each CO or H I velocity field was subdivided into rings 
using the adopted orientation parameters, and a third- 
order Fourier series was fitted to the velocities in each 
ring (cf. Eq.|3|with n=3). The program employs a least- 
squares fitting procedure (singular value decomposition) 
rather than a direct Fourier expansion since the points 
may not be uniformly distributed in azimuth. No extra 
weight is given to points near the major axis, and in fact 
it is points near the minor axis that provide the most 
information about non-circular motions. 

The si and s3 coefficients are plotted against radius for 
each galaxy in Figures|HJa)-^fa), after normalization by 
ci ~ w*, which is also plotted. The error bars represent 
the formal least-squares errors (corrected for the number 
of pixels per beam), using the sum of fractional errors for 
the ratios. While the ci and C3 terms are also affected by 
elliptical or spiral streaming, they are also quite suscep- 
tible to beam smearing effects and, in the case of the ci 
term, dominated by the overall circular rotation. Fluc- 
tuations in the cq, C2, and S2 terms, on the other hand, 
reflect an error in the kinematic center (in which case the 
deviations should fall as R~^) or the presence of an to=1 
mode (lopsidedness) in the potential. Although lopsid- 
edness is clearly revealed in several galaxies, it is usually 
at a fairly low level (|c2| ^ 5 km s"""^), except in the outer 
H I disks of strongly warped galaxies. A more thorough 
analysis of this phenomenon is beyond the scope of this 
paper. 

The lower panels of FigureslHtll2lshow. for each galaxy. 
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Fig. 6. — (a) ci, si and S3 coefficients plotted as functions of radius for NGC 4321. The si and S3 coefficients liave been normalized 
by ci. The CO data is indicated by filled circles and the H I data by open squares. Text labels indicate whether positive values of si 
correspond to radial inflow or outflow. Vertical dashed lines delineate candidate inflow regions given in Table|2 (b) Morphological position 
angle and ellipticity as a function of radius from the isophote fits. In the upper panel, the kinematic minor axis position angle (shifted by 
7r/2) is overplotted as a solid line for the CO and a dashed line for the HI, for all points with estimated error < 5°. In the lower panel, 
the kinematic inclination, expressed as an ellipticity using e = 1 — cosi, is similarly plotted for points with error < 0.05. Corresponding 
values of i in degrees are given on the right axis. The dotted horizontal lines represent the adopted values for the disk position angle and 
inclination, (c) Comparison of si and S3 coefficients normalized by ci, for CO and (d) HI. The heavy long-dashed line represents a linear 
least-squares fit to the data points (with slope given in the upper right corner), while the short-dashed line ("warp line") represents the 
predicted relation between si and S3 for an error in the position angle Fq. The sign of si corresponds to inflow or outflow as indicated at 
the bottom of each plot. 
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Fig. 7.— Same as Fig. [HI but for NGC 4414. 



Fig. 8.— 



Same as Fig. [HI but for NGC 4501. 
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Fig. 9.— Same as Fig. [HI but for NGC 4736. 



Fig. 10.— Same as Fig.lHJ but for NGC 5033. 
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Fig. 11.— Same as Fig.lHl but for NGC 5055. The box in panel 
(d) indicates the region of possible inflow between i?=50"-150". 



Fig. 12.— Same as Fig. [HI but for NGC 5457. 
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the measured si and S3 terms plotted against each other, 
normalized by ci ~ VcS'mi. Coefficients for the CO and 
H I data are plotted separately as panels (c) and (d), 
since they are derived at different spatial resolutions. 
The slope of a linear least-squares fit (heavy long-dashed 
line) is given in the upper right corner of each plot; we 
consider slopes of < 0.1 in absolute value, corresponding 
to Isi/sal > 10, to be good candidates for radial flows. 
The dashed line through (0, 0) is the warp line — the lo- 
cus of points along which an error in the disk position 
angle (F) should fall, assuming a fiat rotation curve. A 
tendency for the data points fall close to this line would 
be consistent with a warp. In addition, a shift in the 
assumed P.A. (Fq) would be equivalent to a shift of all 
data points in a direction parallel to this line. 

If the Si term is interpreted as radial fiow, its sign cor- 
responds to inflow or outflow depending on the sense of 
rotation in the galaxy. Assuming that spiral arms are al- 
ways trailing, we infer that all galaxies besides NGC 4736 
and 5055 rotate counterclockwise in the sky, so that ex- 
cept for these two galaxies the criterion si > corre- 
sponds to radial outflow, as labeled at the bottom of 
each S1-S3 plot. We discuss results for each galaxy indi- 
vidually in t|5.3l 

5.2. Comparison with Isophotal Data 

Figures ICTbY- tT^ b') show, as functions of radius, the 
morphological position angle and inclination, as derived 
from the isophote fits, and the kinematic minor axis posi- 
tion angle and inclination, as derived from the harmonic 
decomposition of the CO and H I velocity fields. Hor- 
izontal dotted lines are used to indicate the adopted 
values of the position angle and inclination from the 
ROTCUR analysis (E^, except for NGC 4736 where 
the adopted inclination is based on the photometry of 
[Mollenhoff. Matthias. & Gerhard (1995). The inclina- 
tion i has been expressed in terms of the ellipticity of the 
projected disk, e = 1 — cosi. Note that the morphological 
inclination, emorph, is expected to decrease in the inner, 
bulge-dominated region of the disk as the isophotes be- 
come more circular. This general behavior is seen most 
clearly in NGC 4414, 4501, 5033, and 5055, which are 
sufficiently inclined so that the fitted values of Fmorph 
and Cmorph are dominated by geometry rather than non- 
axisymmetric structure. 

The kinematic values were derived as follows: 

1. We determined the kinematic minor axis Fmin from 
the harmonic coefficients for each ring by setting 
Vios = Vsys in Eq. 13 and deri ving the roots num eri- 
cally using Ridders' method l|Press et al.ll992l and 
references therein). There are two roots spaced by 

TT, corresponding to the two halves of the axis; we 
adopted the average of the two values (first offset- 
ting one of them by tt) and took half their difference 
as the estimated error a. Only points with cr < 5° 
are shown. 

2. The kinematic inclination or "ellipticity" is defined 
by 

Ekin = 1 - gkin = 1 - g ^1 - (15) 

where q is the cosine of the adopted inclination used 
in the harmonic decomposition. This approxima- 



tion is only accurate to first order in Sq/q, and ad- 
ditional errors can arise if the rings are sufficiently 
tilted with respect to the assumed inclination that 
points in the velocity image get assigned to the 
wrong ring. We checked that our results were con- 
sistent with tilted-ring fits in which the inclination 
is allowed to vary, although such fits generally give 
poorer results because of correlations between fit 
parameters. Only points with a < 0.05 (with errors 
derived from the harmonic coefficients) are shown. 

Based on the results of H4.6I we expect that ellipti- 
cal streaming should be accompanied by a change in the 
morphological position angle Fmorph in a direction op- 
posite to the change in Fmin, assuming the gas orbits 
are aligned with the isophotes. An anti-correlation of 
this type is seen in several of the galaxies, most notably 
NGC 4321, 4501, and 4736, as discussed below. The 
presence of a bar, warp, or strong spiral structure would 
be expected to affect the morphological inclination (via 
the isophote ellipticity) and quite possibly the kinematic 
inclination (via the C3 term) as well. 

5.3. Results for Individual Galaxies 

NGC 4S21 (Ml 00) 

The CO velocity field of this well-studied Virgo cluster 
galaxy shows a strong si term peaking between R—20" 
and 40" with si/ci « 0.35 and a slope \ds3/dsi \ ~ —0.05, 
making it a clear candidate for radial flows. Interpreted 
as radial flow, the si term implies an outflow at ^--^60 
km s""'^ in the plane of the galaxy. However, elliptical 
streaming must be considered as well, since NGC 4321 
is kno wn to possess a stell ar bar with a radial extent of 
~60" l|KnaDen et alJlT995l) . Further evidence for ellipti- 
cal streaming comes from a comparison of the kinematic 
and morphological position angles: the CO kinematics 
show a positive offset in Fmin over the range i?=20"-40" 
[Fig. Elb)], corresponding to a negative offset in Fmorph 
over the same region. Moreover, the H I velocity field 
from i?=40"-180" shows a negative slope in the S1-S3 
plane as expected for a bar potential [Fig. Efd)]. The 
presence of a dominant si term in the CO kinematics 
would be consistent with our modeling of gas flow in the 
ILR region of a barred galaxy [cf. Fig.l^Ja)]. 

NGC UU 

The CO data show nearly circular kinematics, aside 
from a region near R ~ 30" where ds^/dsi ~ —0.1 
[Fig. CJc)], a possible indication of radial inflow. The 
peak inflow speed would be ^8 km s""'^, less than 4% 
of the circular speed. On the other hand, the error- 
bars are large and do not rule out values of ds^/dsi that 
would be characteristic of an elliptical streaming model. 
While there is little evidence for a bar in the i^-band 
image (Fmorph varies by < 2° over the region i?=20"- 
50"), low-amplitude s piral structure had been no ted in 
the if-band image bv lThornlev fc Mundvl l|1997bf) after 
removing an axisymmetric component from the light dis- 
tribution. Thus the apparent inflow signature may be 
related to weak spiral streaming. 

While the H I data show a strong decrease in the si 
component, especially for R > 150", it is well correlated 
with changes in S3 [Fig. [7{d)], suggesting the presence 
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of a strong warp. The slope of the least-squares fit is a 
good match to the warp line, although there is a constant 
offset from that line which may indicate the effects of 
other non-circular motions. A decreasing trend in the 
morphological inclination supports an interpretation in 
which the warp extends as far in as i? w 30", well within 
the optical radius {R25 ~ 100"). This would explain the 
slight mismatch between the CO and H I rotation curves, 
since the average CO and H I inclinations would actually 
be different. We note, however, that effects other than a 
warp, such as finite thickness or flaring of the disk, might 
also affect the morphological inclination. 

NGC 4501 (M88) 

The CO data show an outflow signature very similar to 
that seen in NGC 4321, with implied outflow velocities 
of up to -^45 km s^^ in the region R < 40" [Fig. IH^a)]. 
As in the case of 4321, however, this could well be due to 
elliptical streaming, a possibility that is strongly favored 
by the remarkably clear anti-correlation between Fmin 
and Fmorph in this region. It is unclear whether such 
a signature could be attributed to spiral arms — strong 
spiral features are seen in both CO and Ha maps (Pa- 
per I), but at an inclination of 64° the projection effects 
are severe. The S1-S3 plot for the H I data [Fig. |HId)] 
shows a marginally negative slope at larger radii (with 
si > 0, as would be the case for outflow), which is also 
likely to be a consequence of the bar. Note that al- 
though the isophotes in this outer region (i?=50"-100") 
are nearly circular, isophotes become less sensitive to 
non-axisymmetric structures as the disk inclination in- 
creases. We conclude that NGC 4501 is a likely example 
of a galaxy that is optically classifled as unbarred but 
whose gas kinematics and inner near-infrared isophotes 
strongly suggest the presence of a bar. 

NGC 4736 (M94) 

This prototypical ringed galaxy cont ains a nuclear 
stella r bar and a corresponding CO bar ijWong fc Blitd 
I2000() . We attribute the large Si and S3 terms in the 
CO data to the effects of the bar, although the scatter 
in the S1-S3 plot is large [Fig. I^fc)], partly due to the 
complicating effects of a strong spiral arc at i? '--^ 25" 
whose streaming motions are distinct from the bar's. 
The H I kinematics between R—50" and 150" show a 
strong si term with jds.^/cisi I « 0.2. In an earher paper 
llWong fc Blitd l200(l . we had argued that the relative 
weakness of the S3 term appeared to rule out an ellip- 
tical streaming model. The more sophisticated analysis 
presented in H4.2l suggests that instead we may be view- 
ing elliptical streaming in the vicinity of an ILR; in fact, 
the prominent Ha ring in this galaxy at i? = 45" has 
been modeled as the ILR of a l arge-scale oval distortion 
l|Gerin. Casoli. fc Combes 1991) . A pronounced decrease 
in Fmorph occurs between i?=50"-100", also consistent 
with such a distortion [Fig. El^b)] . While a warp model 
is not ruled out by the kinematics alone, due to the low 
inclination of the galaxy, the observed anti-correlation 
between Fmorph and Fmin would then be harder to ac- 
count for. 

NGC 5033 



The CO data show some evidence for a warp based 
on the S1-S3 plot [Fig. EIc)], although the uncertain- 
ties are large. No indication of radial inflow at speeds 
of >5% of the circular speed (wc ^ 200 km s~^) is seen 
for R < 55". For the H I data, the si and S3 terms 
are positively correlated at radii beyond 90", raising the 
possibility that the galaxy's warp, which only becomes 
prominent around i?25=320", is pervasive throughout 
the galaxy, as suggested for NGC 4414. Such a distur- 
bance may be related to the weak interaction of NGC 
5033 with Holmberg VIII, as evidenced by the tidal H I 
plume on the southern side of the galaxy l|Thean et alJ 
1199 7^. However, the peak in si and S3 near i?=110" 
[Fig. ITUT a')] might instead be the result of spiral stream- 
ing. The morphological inclination also falls with radius 
beyond i?=30"-100", although it does not closely follow 
the kinematic inclination [Fig. llOf b)] . It is possible that 
extinction or projection effects due to the galaxy's high 
inclination may dominate over the warp, as far as the 
isophotes are concerned. A deeper K-hand image would 
be useful for clarifying the situation. 

NGC 5055 

A strong negative slope in the S1-S3 plot for the 
CO data [Fig. [TTT c)] suggests the presence of low-level 
streaming motions centered around _R=30". Although 
inspection of the i^'-band image does not reveal clear ev- 
idence of a bar, the isophotes do become "boxier" around 
this radius, leading to a dip in the ellipticity Cmorph- In 
spite of this localized deviation, the smallness of the si 
term places an upper limit of vr < 5 km s~^ (3% of the 
circular speed) on the magnitude of any radial flows in 
the region i?=15"-70". 

The H I kinematics are clearly dominated by a warp 
for R > 200" (0.5i?25)- Inside of this radius, a possi- 
ble region of gas inflow is seen, highlighted by the small 
box in Fig. Ilir dV We interpret this weak feature (again 
amounting to radial motions of no more than 3% of the 
circular speed) to spiral streaming motions, since one 
finds CO and H I spiral arms crossing the galaxy's minor 
axis at approximately these radii, where their effects are 
more likely to have an impact on the observed kinemat- 
ics. Fluctuations in the isophotal parameters Fmorph and 
emorph are also apparent around R ~ 70" [Fig. [njb)], 
and are likely associated with spiral arms. 

Note that for this galaxy, emorph never rises to a value 
comparable to ekin- Aside from a few humps that may be 
attributed to spiral structure, the morphological inclina- 
tion reaches a fairly constant value of around 57°, while 
the kinematics favor an inclination of 63°, also with very 
little scatter (this comparison is made well inside the re- 
gion of the warp). One possibility is that the morpholog- 
ical inclination is reduced by the finite thickness of the 
stellar disk. A radius of 200" corresponds to ^10 kpc 
for this galaxy, so the projected minor axis at this radius 
should have a length of ~4.5 kpc if the true inclination is 
63°. If the minor axis is extended by 1 kpc due to a thick 
disk, however, the morphological inclination is reduced 
to 56°, consistent with the observed value. 

NGC 5457 (MlOl) 

The inclination of this galaxy is not well determined 
by either morphology or kinematics; the adopted value 
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Table 4 

Candidate Regions For Radial Flows, Based On 

Sl/S3 



Galaxy 


Tracer 


Radial Range 


Maximum vn 


Interpretation 


NGC 4321 


CO 


20" 


-40" 


56 ± 5 


Bar streaming 


NGC 4414 


CO 


15"- 


-45" 


-8± 2=^ 


Spiral streaming 


NGC 4501 


CO 


15"- 


-40" 


48 ± 5 


Bar streaming 


NGC 4501 


HI 


40"- 


150" 


10 ±4 


Bar streaming 


NGC 5055 


HI 


50"- 


150" 


-6± 1 


Spiral streaming 


NGC 5457 


HI 


250"- 


-600" 


-56 ±5 


Outer warp 



"Negative values denote inflow. 



is based on fitting just two rings in the H I velocity field. 
However, our value of 21° is close to the 18° value of 
RC3 as we ll as the 22° tha t resul ts from the Tully-Fisher 
relation of iPierce" fc Tullvl l)1992l) using an H I lincwidth 
of 190 km s-^ 

Despite poor signal-to-noise in the inner arcminute, the 
CO velocity field shows strong evidence for large si and 
S3 terms [Fig. I12f c')]. and indeed a substantial twist in 
the kinematic minor axis is seen in the tapered CO veloc- 
ity field (Fig. The slope of the S1-S3 relation appears 
to be negative and thus indicative of elliptical streaming. 
Identifying the bar in the optical images is difficult due to 
the strong spiral structure in this galaxy, but the bar-like 
CO and Ha morphology, as well as the rise in the ellip- 
ticity of the isophotes across the region _R=20"~-60" to a 
value much larger than in the rest of the disk, confirm 
the presence of a weak bar or oval distortion in the inner 
regions. The H I kinematics show a strongly decreas- 
ing si term beyond i?=200" (0.2i?25) which is probably 
dominated by a large outer warp, with additional fluc- 
tuations in si and S3 at certain radii due to the strong 
spiral structure in this galaxy. Unfortunately, the incli- 
nation is too low to clearly distinguish the effects of the 
warp from possible radial flows. 

6. DISCUSSION AND CONCLUSIONS 

We have carried out a search for radial flow signa- 
tures in the CO and H I velocity fields of seven nearby 
spirals by decomposing each elliptical ring into a third- 
order Fourier series. By requiring that jsi/saj > 10 over 
a significant range in radius, where si amd S3 are the 
coefficients of the smip and s'mSip terms, we identified 
candidate regions for radial inflow or outflow as summa- 
rized in Table 01 Aside from NGC 5457, where an inflow 
signature is degenerate with a warp, we flnd photometric 
evidence for bars or spiral structure in all such regions, 
suggesting that elliptical streaming in a bar or spiral po- 
tential is the dominant contributor to non-circular mo- 
tions. While inflow may be superposed on these motions, 
we find no unequivocal evidence for radial inflow alone. 
Three of the galaxies, NGC 4414, 5033 and 5055, show 
nearly pure circular rotation in their inner few kilopar- 
secs (i? < 60"), with an upper limit to any radial inflows 
of ~5-10 km s~^, or 3-5% of the circular speed. These 
are probably the strictest limits that can be placed on 
radial gas flows in external galaxies. 

Although our radial inflow limits are a small fraction 
of the circular speed, they cannot be used to rule out 



theoretical models that invoke radial flows, since such 
models generally predict even smaller inflow speed s. On 
the higher end of estimates is that bv iBlitd l)1996f) . who 
estimates that an inflow velocity of —7 km s~^ would be 
needed to completely replenish gas consumed by star for- 
mation in the inner Galaxy with H I from twice the solar 
radius. Given that estimates of the gas consumption time 
in our sample of galaxies are substantially longer than in 
the Milky Way (P aper I), such a high inflow rate would 
not be expected. iStruck-Marcelll (|1991l) discusses how 
radial flows can serve to maintain hydrodynamic stabil- 
ity in gaseous disks, but such flows are always subsonic 
(<5 km s~^) and in fac t reduce to ze ro for an gas 
surface density profile. iStruck &: Sm ith (1999) further 
develop this idea into a steady-state model for a turblent 
multiphase ISM, predicting slow (~3 km s^^) inflows of 
cold gas in the galaxy midplane. Pure chemical evolu- 
tion mode ls invoke even smaller infl ow speeds (e.g. ~1 
km s~^ in lPortinari fc Chiosil 1200(11) since they are con- 
cerned with steepening metallicity gradients in galaxies 
over timescales of many Gyr. Thus, our observations and 
analysis techniques have yet to reach the regime where 
such flows can be easily detected. 

The results presented in this paper underscore the dif- 
ficulty of detecting radial inflow in disks where a warp or 
m=2 distortion in the potential exists. Once such effects 
become dominant, they cannot be easily fitted and re- 
moved without performing a more sophisticated analysis 
than that performed here. For instance, by neglecting 
the important role of shocks in the gas component, we 
are unable to model the very processes (bars and spiral 
arms) that are most likely to lead to radial inflows in the 
inner regions of galaxies. Even with more sophisticated 
modeling, however, a leap in sensitivity and angular res- 
olution will be needed to examine shock regions in detail. 
Similarly, our lack of knowledge about warps limits our 
ability to analyze the gas kinematics in the outskirts of 
galaxies. The default assumption, that warped orbits re- 
main circular even as they move away from the principal 
plane, must be only a crude approximation. Since warps 
may be intimately related to infalling gas (Binncy 1992), 
which in turn could induce radial flows, further modeling 
is clearly desirable. 

An intriguing result from our analysis has been the 
possible detection of low- level kinematic warp signatures 
in NGC 4414 and 5033, at radii well inside the radius at 
which the warp becomes prominent (which is roughly the 
optical radius -R25)- While spiral structure may confuse 
the picture in NGC 5033, it is unlikely to play a major 
role in NGC 4414, where the spiral structure is much 
more flocculent and the isophote fits seem to indicate an 
optical warp as well. The prominent kinematic warps in 
NGC 5055 and 5457 also app ear wel l inside R25, contrary 
to usual expectations (e.g., lBriggal [i990). Since warps 
are often straightforward to identify in the S1-S3 plane, 
the harmonic decomposition technique will be valuable 
for assessing the radius at which warping begins, which 
in turn may shed light on the physical mechanisms that 
maintain warps. 
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APPENDIX 



EFFECT OF DISSIPATION ON GASEOUS ORBITS 

For a non-zero pattern speed, a coUisionless weak bar model has the property that the stable period ic orb its change 
orientation (from parallel to perpendicular to the bar or vice versa) with each resonance crossing [Figure UTsT a')]. Many 
of the resulting orbits intersect each other, and are hence not accessible to the gas component, since only one velocity 
is permitted for a fluid at a given spatial location. Instead, the gaseous orbits are expected to turn gradually between 
resonances, producing a spiral pattern [see Figure IXTsT b)]. Since cloud collisions, shocks, and energy dissipation lead 
to vastly more complex behavior than can be described by a model such as the one given in ij4.2l this situation is 
generally examined with the help of numerical simulations, employing either a hydrodynamical or "sticky particle" 
approach. 

One can obtain an approximate analytic description of gaseous orbits in a b arred potentia l by intr oducing a damping 
term into the equ ations of motion, as has been done by various authors (e.g.. iLindblad & L indblac]l ll994tlWadalll994 
iBvrd et al.lll998|) . The additional term simulates a frictional force that is proportional to the velocity perturba tion in 
the ra dial and/or azimuthal direction. For simplicity we consider damping in the radial direction onl y, following I Wadal 
l(T99l . (The more inv olved case of both radial and azimuthal damping is treated bv iBakeil 1200(1) . The first-order 
equation of motion (cf. iBinnev fc Tremaiii3ll987l p. 148) is then: 



Ri + 2Ai?i + kIRi 



dR R{n - np) 



COS mfpo 



(Al) 



-Ro 



where A is the damping t erm a nd = (^o ~ ^p)t. The solutions for (i?, </>, w^, w^) in the rest frame of the bar are 
given bv I Sakamoto et all l)1999|) . although we have redefined some of the constants here to facilitate comparison with 



R — i?( 



1 - — cos(2(/)o + oq] 
B 



00 + J sin(2(?!)o + /3o) 
VR^Vc{l - ujp)[A sin(2(?!)o + 5q)] 

V^=Vc{l - LOp) 



A B 
1 - - cos(20o + ^o) + y cos(2(/)o + /3o) 



The amplitudes of the epicyclic motion are given by 

2 



A- 



2$2(i?o) 



B^^/{E + FY + 2EF{cos So - 1 



where 



E 



2A 

1 — LJ„ 



The phase shifts of the orbit are given by 

tan So — - 



4A(nn - 



tan f3o 



i?o(l - i^p) 
2^2 (i?o) 

vui~u;py ■ 

sin Jo 



$^(i?o) 



(A2) 

(A3) 
(A4) 

(A5) 

(A6) 
(A7) 



(A8) 



Ao cosSo + iF/E) 

Note that in the limit that A ^ 0, we have A ^ a2 and B (02 + 62), in terms of the coefficients of the coUisionless 
model. 

Transforming to the galaxy rest frame, we find that for t = 0: 

VR^Vc[A{l-tJp)sm{2e + So)] (A9) 



ve^Vc 



A B 
1 ~ - cos(26' + So) + —(1 - LOp) cos(26i + Po) 



(AlO) 




Fig. A13. — Orbits viewed in the rotating frame for a (a) dissipationlcss bar model and (b) dissipative bar model with a flat rotation 
curve. The dashed lines at radii of 1 - v^/2, 1, and 1 + ^2/2 are the ILR, CR, and OLR respectively. 



Recall that the harmonic coefRcients are defined by: 

V\os = Vsys + {vg cos V' + v_R sin V') sin i 

— cq + ci cosip + Si sin-0 + cs cos 3-0 + S3 sin3'0 , (All) 
where ip — 9 — 0obs — 7r/2. Thus we find: 

CO^Vsys (A12) 

— = 1 - ][{-'2ujp + l)v4cos(26iobs + So) + (1 - t^p)Bcos(26iobs + f3o)] 

— = i [{-2ujp + 1)A sin(20obs + So) + (1 - u;p)B sin(2^obs + M] 

— = -i[(2Ljp - 3)Acos(2^obs + 60) + {I- Wp)-Bcos(20obs + Po)] 

— = j[(2wp - 3)Asin(20obs + ^o) + (1 - c^p)Ssin(20obs + f3o)] 

where = VcSini. These relations give the observed harmonic coefficients in terms of the orbit parameters 
iA,B,(3o,So). 

For the case of an elongated potential discussed in M.2\ the orbit parameters {A, B, 60, Po) are given by: 

A^ 'P°* — (A13) 

(1 - u;p)^{l ~ 4iOp + 2a;2)2 + 4^2(1 „ ^^)2 

Epot V 4Aepot 



B = -\ T^ + TT^ +jf^{cosSo~l) (A14) 



2A(1 - Wp) 
1 - 4wp + 2ujI 



tan^o^ /-r."T.2 (A15) 



tan/3o = '^^".f" . TT . (A16) 

cos(5o + epot/[2A(l-cjp)] 

where A = A/fio- With our choice of signs, So — tt at ujp=0 for A = 0. Model orbits calculated using these equations, 
using a flat rotation curve, are shown in Figure lA13r b'l. Note that the damping term leads to a gradual change in 
orbit orientation near the ILR {uip = 1 — \/2/2), although the orbit oscillations still go nonlinear near the CR {uip = 1). 
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